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Abstract 

We present a formulation of the operator product expansion that 
is infrared finite to all orders in the attendant massless non-Abelian 
gauge theory coupling constant, which we will oftentimes associate 
with the QCD theory, the theory that we actually have as our primary 
objective in view of the operation of the LHC at CERN. We make 
contact in this way with the recently introduced IR-improved DGLAP- 
CS theory and point-out phenomenological implications accordingly, 
with an eye toward the precision QCD theory for LHC physics. 

1 Introduction 

With the start-up of the LHC the era of precision QCD, by which we mean 
predictions for QCD processes at the total precision tag of 1% or better, 
is upon us and the need for exact, amplitude-based resummation of large 
higher order effects is becoming more and more acute. Methods to facilitate 
the realization of such resummation are then of particular interest. In this 
paper, we revisit the pioneering use of operator product expansion (OPE) 
methods, as presented by Wilson [lj for short- distance limits of physical 
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processes and as applied by Gross, Wilczek and Politzer in the QCD j2] 
theory, especially as it is realized in the DGLAP-CS [2111] theory, from the 
standpoint of resummation of its large infrared effects with an eye toward 
the attendant application of the corresponding parton model representation 
to LHC precision physics. In this way, we make contact as well with the 
recently introduced IR- improved DGLAP-CS theory in Refs. [5HZ]. 

Specifically, it is well-known [HHID] that the usual formulation of the Wil- 
son expansion in massless gauge theory is infrared divergent: the easiest 
way to realize this is to note that, already at one-loop, the respective lead- 
ing twist operator matrix elements between fundamental particle states are 
in general infrared divergent and must be evaluated at off-shell (Euclidean) 
points in massless gauge theory - see for example Refs. [5HTU]. The result 
is that the coefficient functions of the Wilson operators in the OPE which 
encode the leading Q 2 dependence of the expansion are in general infrared 
divergent order-by-order in renormalized perturbation theory^. Of course, all 
such infrared divergences cancel in physically observable (hadronic) matrix 
elements of the expansion so that, from the standpoint of such observables, 
the issue is one of choosing the best rearrangement of the large infrared effects 
that remain after all infrared divergences have canceled. Here, we will resum 
these large infrared effects. As a result, in what follows, we reformulate the 
OPE in such away that the respective expansion components are infrared 
finite. As a further result, we show how the new IR-improved DGLAP-CS 
theory in Ref. [5H7] arises naturally in this context. We argue that the IR- 
improved expansion should be closer to experiment for a given exact order 
in the loop expansion for the coefficient functions and respective operator 
matrix elements. 

The paper is organized as follows. In the next section, we recapitulate 
the formulation of the OPE following the arguments of Wilson as used in 
Refs. [&MTU] for the analysis of deep- inelastic lepton-nucleon scattering [TT] . 
the proto-typical physical application of the method. In Section 3, we show 
how to improve it so that its hard coefficient functions are IR finite. We also 

1 In Ref. pp, Wilson pointed-out already that the coefficient functions in his expansion 
could be calculated order by order in perturbation theory and that the n-th order term 
could contain logarithms of z 2 m 2 where z is the space-time interval of the respective two 
operators and m is the free field mass, so that these logs would be divergent at m = 0. 
He also noted that an arbitrary subtraction constant a could be introduced to convert the 
argument of these logarithms to z 2 a 2 . This is equivalent to what we have stated in the 
text with a = [i where /i is identified as a Euclidean point in an appropriate convention. 
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make contact with the new IR- improved DGLAP-CS theory [SH7] . In Section 
3, we also sum up with an eye toward phenomenological implications. 



2 Review of the OPE 



For pedagogical reasons, we follow the historical development and use the 
deep inelastic electron-proton scattering problem discussed so effectively by 
Bjorken [12J as our starting point: e~ (£)+p(p p ) — > e~ (£')+X(px)- Indeed, his 
discussion set the framework for the issues we address here. The kinematics 
and notation are summarized in the Fig. [H so that we use x = Xsj = 
Q 2 1 \2m p v) for Bjorken's scaling variable which has the interpretation in the 
attendant parton model as the struck parton's momentum fraction when 
v = qpp/mp with q = £ — £', Q 2 = —q 2 - In the Fig. [TJ the parton momenta 
are Pi(p'i) before(after) the hard interaction process. The limit of Bjorken is 



Figure 1: Deep inelastic electron-proton scattering: q = £ — £', v = 
qPp/rrip, x = XBj = ~q 2 /(^ m P ^), £(£') is the four-momentum of the ini- 
tial(final) e~, pa is the four-momentum of A, A = a,p, where a is a parton. 

then of interest here, in which we take Q 2 — > oo with x fixed. In this limit, 
where here we will for reasons of pedagogy focus on the photon exchange in 
Fig. [DIE the standard methods can be used to represent the imaginary part 

2 As it is well-known, adding in the effects of the Z exchange is straightforward and 
does not require any essentially new methods that are not already exhibited by what we 
do for the photon exchange case. 




P(Pp) 
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of the attendant current-proton forward scattering amplitude as 

W*e M ( Pp ,q) = ± I d'ye^<p\[jf M (y),J EM (0)}\p> 

(-gap + qaqp/q^W^q 2 ) (1) 

+ —(p p - qqp P /q 2 )a{pp - qqp P /q 2 )pW 2 {v, q 2 ) 

m p 

Here, </<f M (?/) is the hadronic electromagnetic current and Wip are the usual 
deep inelastic the structure functions, which first were shown to exhibit 
Bjorken scaling by the SLAC-MIT experiments [11] already at Q 2 = l + GeV 2 , 
precocious scaling - we return to this point below. For our purposes here, 
henceforward we drop the superscript on J EM so that J EAI = J for ease of 
notation and we always understand the average over the spin of the proton 
even when we do not indicate so explicitly. In Bjorken's limit, we have 

lim m p W\ (u, q 2 ) = Fi(x) 

Bj (2) 
limz/H/ 2 (z/,g 2 ) = F 2 (x) K ' 

Bj 

where the scaling limits F 12 only depend on Bjorken's variable x and we 
denote 

lim = lim \ x ^ &xed . 

Bj 

The QCD theory of Gross, Wilczek and Politzer (2] provides a quantum field 
theoretic explanation of the observed Bjorken scaling behavior via Wilson's 
OPE. 

Specifically, in Bjorken's limit, the phase in integral over space-time in 
(pQ) oscillates rapidly except in regions where it is bounded so that the value 
of the integral is dominated by the latter regions, which are well-known to 
correspond to the tip of the light-cone |13j . the short-distance regime. Using 
Wilson's expansion in this regime, we get the OPE [81410] 

-. oo 

+ :£ £ ~ ^H n (0)y^ ■ ■ ■ + • • • , 

y ley °n=0 j 

(3) 
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where we have neglected gradient terms without loss of content for our pur- 
poses here and as usual e 4- 0. We also note that {Qi (y)} are traceless, 
symmetric spin n operators of dimension n + 2 or of twist = dimension -spin 
= 2 [13J. The ■ •• represent operators of higher twist that are suppressed 
by powers of q 2 to any finite order in perturbation theory. The coefficient 
c-number functions {Cj^ } are dimensionless and can be computed in renor- 
malized perturbation theory. 

Continuing our recapitulation of the methods in Refs. [STfTQ] , if we define 
the spin averaged proton matrix elements of the operators O- 7 via 

< P\Oi v .., n (0)\p > | spin averaged = i n —Pp m ■ ■ -P Pfln M? + -.., (4) 

I iip 

where the second • • ■ denotes trace-terms, we get the following relation- 
ship [8 HIU1H1] between the moments of the structure functions and the Fourier 
transforms of the coefficient functions: 



f dxx n F 1 (x, q 2 ) = V Cf+ l \q 2 )M\ 
dxx n F 2 {x,q 2 ) = Y,C [ ${q 2 )M : 



(5) 



n+2 
3 ' 



where [8] 



2 \ dq 2 J J " y 2 — iey 

The q 2 dependence of the C^ 1 ' is controlled by the Callan-Symanzik equa- 
tion jl] which reads 



SjSV) = W) n+1 (~T [ dSe^^l. (6) 



C# = (7) 



where \x denotes the renormalization scale, 

for the attendant renormalized coupling g, and the anomalous dimension 
matrix ^\g) is given as 

= (^° 1 ^]^ l ^ C ^j Iso.regularization fixed (8) 
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where the operators 0\ n = 



are renormalized via 



(9) 
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so that they mix under renormalization in the well-known way [8 , 9j and we 
use a standard notation of the renormalized, R, and unrenormalized , bare, 
operator representatives, go is the bare coupling constant. It is well-known 
that the solution of ([7]) leads to the conclusion that the asymptotic Bjorken 
limit is controlled by the operators with the smallest eigenvalue for their 
anomalous dimension matrix in an asymptotically free theory such as the 
QCD [2] which we have in mind here. The implied behavior for the RHS 
of (JSJ) is in agreement with experiment Here, we want to focus on the 
IR-improvement of the C-^, M". 



The isolation of the infrared aspects of the is immediate if we use the 
fundamental particles in the respective Lagrangian quantum field theory, 
quarks and gluons in the case of QCD, to evaluate the essential anomalous 
dimension matrix elements J^ig), as this is equivalent to studying deep in- 
elastic scattering from these fundamental particles and takes us immediately, 
at least conceptually, to the parton model perspective studied famously by 
many [BH1I]. 

Specifically, we then focus on the parton level version of hadronic tensor 
W a p which for definiteness we associate with a fermion F in the underlying 
asymptotically free theory(QCD): 



where we use the fact that q° > to drop the remaining term in the com- 
mutator and we always average over the spin of the fermion F, as we do for 
the protonp. We see clearly from ([T]) that the RHS of ffTOl and that of (pP) 
involve the same OPE. 



3 IR-Improved OPE 




— / dV® <PF\[Jp(v),Ja(0)]\PF> 

{2Txf^5(q + p F -p x ) <Pf\M0)\px >< Px\J a {0)\p F > 



x 



(10) 
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We first focus on the matrix element 

M X ,a =< Px\Ja(0)\p F > • (H) 

Following Refs. [SIE], we isolate the dominant virtual IR divergences associ- 
ated to the incoming line via the formula 

M X ,a = e°° ^ CD < Px\Ja(0)\PF >IRI-virt, (12) 

where the virtual infrared function Bqcd is given in Refs. [5j[6] . The RHS of 
this last equation is valid to all orders in a s = g 2 /(Air) so that one computes 
< Px\J a (Q)\PF >mi~virt from < Px\J a (®)\PF > by comparing the coefficients 
of the powers of a s on both sides of (j!2p iteratively. 
Introducing this result into ffTUl) . we arrive at 

W£p( PF ,q) = (2^) 3 ^% + PF - p x )e 2a ^ B ^ D iRi- vir t< Pf\M0)\ P x > 

x (13) 

< Px\Ja{0)\pF >IRI-virt ■ 

We next isolate the leading soft, spin independent real emission infrared 
function associated to the incoming line as follows. We first separate {X} 
into its multiple gluon subspaces via 

{X} = {X : X = X' <g> {Gi <g) . . . <g) G n }, for some n > 0, X' is non-gluonic}. 

(14) 

Then we have 



e 2as ^ CD IRI- Vir t< Pf\M0)\ PX > < Px\Ja(0)\p F >IRI-vi 



virt 



a 2a s ^B QCD 



SQCD(kl) ■ ■ ■ Sq C d(K) IRI-virt < Pf\J p{^)\PX' > 

< PX'\Ja{0)\pF > IRI-virt H h IRI-virthreal < PF \ J/3 (0) \px' , h , ■•• , k n > 

< px',h, ■ ■ ■ ,k n \J a (0)\p F >IRI —virt&zreal j 

(15) 

where the real infrared function S>QCD(k) is given in Refs. [SJE]- The IR- 
improved quantities 

I R.I —virtlkreal 

< Pf\Jp{0)\px > < Px\J a {0)\pF >IRI —virtSzreal 

are defined iteratively from f|T2l) . f|T5|) to all orders in a s and they no longer 
contain the infrared singularities from the initial line associated to Bqcd and 
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to Sqcd, although, because of the non-Abelian infrared algebra of the theory, 
they do contain other IR singularities which of course cancel in the structure 
functions by the KNL theorem for massless fundamental fermions. For mas- 
sive fundamental fermions, these latter singularities also cancel provided we 
resum the theory as we are doing here accordingly - see Refs. [21] . 
Introducing the representation in (|15p into ffTUj) we get 

Wjf a (p F , q) = (2tt) 3 % +PF- Px)e 2aMB ^ CD [s QCD (h) ■ ■ ■ S QCD (k n ) 

X 

IRI-virt < PF\Jp(0)\pX> > < PX'\Ja(0)\pF >IRI-virt H 

+ IRI-virtkreal < PF\ Jp(v) \P X', k\, • ■ ■ , k n > 

<Px',h,--- , k n \J a (o)\p F > IRI 

—virt&ireal 

2vr J y ^^n\J 1=1 k° 

J X' n J 

I RI —virt&real 

< Pf\J/3{0)\px', h, ■ ■ ■ , k n > 

< PX>, h,--- , k n \Ja{0)\pF >IRI —virt&ireal 



1 

2^ 



rf4y e iqy e SUM IR (QCD)+D QC D 
I RI —virtlkreal <PF\[Jp{y),Ja{0)]\PF > I RI —virt&creal t 



(16) 



where we have defined 

SUM IR (QCD) = 2aMB QC D + 2a s B QCD (K m&x ) 

d 3 k~ 



/d k ~ 
-r^-S QC D(k)8(K max - k), 

/d 3 k ~ 
—S QC D(k) [e"** - 9(K max - k)] , (17) 

and we stress that (fTBl) does not depend on K max . Using the standard par- 
tonic view, by which we have 

W ^ = E/ (18) 

„ Jo x 
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for appropriately defined parton distribution functions {J- a }, we introduce 
the OPE in ([3]) into (ITB|) and use ([T]) to get the IR-improved results 



f dxx n F 1 (x, q 2 ) = V ^ +1) (g 2 )M; +1 
Jo j 

[ 1 dxx n F 2 (x,q 2 ) = J2^ n 2(<l 2 )M, 
Jo 



(19) 



n+2 
3 ' 



where [S] 



Cm<?) = -i(q 2 ) n+ ~tto / d 4 ye iqy e bUM ^ UU)+u ^ D ^ w ' (20) 

J ' 2 \ dq 2 J J y 2 -im 



and now 



Pl^ui—u n (0)|P ^ Ispin averaged 



IRI-virtkreal < P\Oj il ... /in (0)\p > IRI-virtkreal | spin averaged — Pp w ' ' 'Pp^ n ^j 



(21) 

where the second • • • again denotes trace-terms and the {C^} are the re- 
spective (new) IR-improved OPE coefficient functions. The q 2 dependence 
of the is also controlled by the Callan-Symanzik equation |3] which now 
reads 



C«"> = (22) 



where now the new matrix ^\g) is determined by the renormalization prop- 



erties of the IR-improved matrix elements in (}2~T]) as we will discuss presently. 
We need to stress that in writing ( 120]) we work to one-loop order in the various 
coefficients in this paper. 

A convenient starting point for obtaining the new matrix j^\g) is pre- 
sented by the pioneering analysis of the authors in Ref. [T21E0] - Working 
directly from the the representation in (fl8l) . the authors in Ref. [19] make 
contact with the matrix 7^ (g) for the unimproved OPE as follows. Focusing 
for definiteness for the moment on the non-singlet operator [8] N F ' b (y) = 
\i N - l Si)(y)-i^V M2 ■ ■ ■ V m \ b ip(y) - trace terms, where V M = + igr a A^ is 
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the covariant derivative, X b is a flavor group generator and S denotes sym- 
metrization with respect to the indices [i% ■ ■ ■ fi n , we have the matrix element 
between fundamental fermion states, where spin averaging is understood here 
as well, as 

< p\ N F >\y)\p >= F '*0 N (.a a , e) Pfil ■ ■ - PflN - trace terms (23) 

where we use d = 4 — e dimensions for regularization and the notation 
F,b O N (a s , e) = Mp b to make immediate contact with the arguments in Ref. [19] 
The renormalized matrix element F,b O N (a s , e) is related to the bare one as 
we have indicated in ([9]): 

F ><V N (a s ,e,p 2 /^ 2 ) = Z Q \a S} ~) F 'toL((a s )(/i 2 /p 2 ) E ,e) (24) 

so that collinear divergences are regularized by taking p 2 ^ in the approach 
in Ref. [19]. Using an arbitrary vector A with A 2 = we get 



F V(a s ,e) =<p\ N F > b .. m (y)\p>A^...A^/(Ap) N . (25) 



The pole part of F ' b O N which is the renormalization part Zq 1 can be de- 
termined in any gauge by gauge invariance. We set A = n where Xsj = 
np/npp, x = kn/np and nA a = so that we are in a light-like gauge. This 
allows us to write, following Ref. [T9] . 



F ' f 7V (« s , e) = J' dxx N ~ 1 F 't>(x, a., e) (26) 



where 



V(x, a s , e) = Z F [5(x - 1) + xj— -/{x - — ) [^-T(p, k) p] (27) 

{2-K) a pn Akn 

where we use the notation of Ref. [19J so that T(p, k) is the respective fully 
connected four-point function and \fi~B denotes ^2 aa , b aa iB^, with corre- 
sponding notation for B $\. Zp is the fermion field renormalization constant 
as usual. By first analytically continuing the LHS of f l24|) to d = 4 + e dimen- 
sions with e > the authors in Ref. [19] note that the limit p 2 — > gives the 
RHS as just Zp 1 (a s , e) for an appropriate normalization of X b . The RHS of 
( |26|) may then related to the moments of the densities of partons in a quark, 
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Ts(x, a s , 1/e) in the notation of Ref. [TU], by writing a dispersion relation for 
[(fl/ (4kn))T(p, k) tf\ and performing the attendant k 2 integral by closing the 
contour around the dispersive poles(see Sect. 4.2 of Ref. [E]), analytically 
continuing to d = 4 + e dimensions with again e > and finally taking the 
limit p 2 — > to get 

If 1 1 1 

Zo 1 (a s ,-)= / dxx N ~ 1 [T qq (x,a s ,-)9(x) -T qq (-x,a s ,-)9(-x)} (28) 

where r^r^) is the respective parton density for a quark ( ant i- quark) in a 
quark. The coefficients of - on both sides of this last equation then give the 
fundamental result, derived in Ref. [T§] . 

- 7 (7V) (a s ) = 2 J dxx N - l [P qq {x,a s )9{x) - P qq -{-x,a s )9(-x)} 
= 2[P qq (N,a s ) + {-l) N P qq (N,a s )} 

where we define 

F(N) = I dxx N - 1 F(x) 
Jo 

and the Pba are the usual DGLAP-CS [Slllj splitting kernels defined in the 
convention of Ref. p3] and ^ N \a s ) is the respective anomalous dimension 
of the operator N F,b . 

To apply this calculation to our new anomalous dimension matrix we IR- 
improve it at each step as we have shown above (and as we have shown for 
the IR- improved DGLAP-CS theory in Refs. [HE]), so that we replace 

< p\ N F ' b (y)\p >^< p\ N F ' b (y)\p > 

as defined in ff2T]) with the corresponding substitution of F,b O N {a s , e) by the 
analogous F,f O N (a s ,e). This leads to the relationship 

F ' b d N (a s ,e jP 2 /f, 2 ) = Zz\a s , \) F ^ are ((a s )^ 2 /p 2 )\ e) (30) 

between the renormalized and bare IR-improved matrix elements. The analoga 
of J25D and (J2ZD are then 

F ' b O N (a s , e) = j 1 dxx"' 1 ^(x, a s , e) (31) 
11 



where 



F > b 0(x, a s , e) = Z F [5(x - 1) + x^S(x - —) [Jr f ^ k ) ( 32 ) 

(2ir) a pn 4kn 

and we continue to use the notation of Ref. [19] so that T(p, k) is the re- 
spective IR-improved fully connected four-point function obtained from the 
unimproved one, T(p, k), by using the master formula Eq.(l) in Refs. [6] re- 
stricted to its QCD aspect, for example. This means that we get the analog 
of (J28D as 

Z 6 \a s ,-) = J dxx N - l [V e ^{x,a s ^)B{x)-V^{-x,a s ^e{-x)\ (33) 

where T^ p , Y e q x ^ are the respective IR-improved parton densities. We get in 
this way the identification of the respective IR-improved anomalous dimen- 
sion as 

- 7 W («s) = 2^[P ? 7(iV, a.) + (-l)^(^ a.)] (34) 

where the Pgq P , Pqq V are the respective IR-improved kernels as introduced in 
Refs. [51E], where we advise that the notation of Ref. [19J differs from that in 
Refs. [HIE] by whether or not one includes the factor a s /(27r) on the RHS of 
( 130]) in the definition of the kernels. This allows us to write at IR-improved 
one-loop level the identifications 

-7 (JV) K)^2g^(AT) (35) 

where the labels i,j span the usual values for the one-loop anomalous di- 
mension matrix for the evolution of the parton distributions as given in 
Refs. [3tmi8l4T0] for example. This establishes in a rigorous way the con- 
nection between the IR-improved DGLAP-CS theory in Ref. [51E] and the 
OPE methods of Wilson as used by Refs. [814T0] in the study of deep inelastic 
lepton-nucleon scattering. 

Evidently, this connection may be manifested in the analysis of other 
physical processes as well. We refer the reader to Refs. [6|,[7] wherein the new 
precision-baseline MC Herwiril.031 which realizes the IR-improved DGLAP- 
CS kernels has been introduced and compared to the Tevatron data [22|[23] 
on single Z production. Its application to the various physical processes at 
LHC is in progress and will appear accordingly elsewhere [24], where we 
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need to stress that Herwiri 1.031 can be applied to any process to which 
Herwig6.5 [25] can be applied and that it interfaces to MC@NLO [26J the 
same way that does Herwig6.5. As we have shown in Refs. |3|7], we have an 
improved agreement between the IR-improved MC's shower and the Tevatron 
data with no need of an abnormally large intrinsic transverse momentum 
parameter, PTRMS~ 2GeV in the notation of Herwig [22], as it is required 
for similar agreement with Herwig6.5 [27] • We point-out that, consistent with 
the precociousness of Bjorken scaling, the IR-improved MC Herwiri. 031 gives 
us a paradigm for reaching a precision QCD MC description of the LHC data, 
on an event-by-event basis with realistic hadronization from the Herwig6.5 
environment, that does not involve an ad hoc hard scale parameter, where 
we define "hard" relative to the observed precociousness of Bjorken scaling. 
What we have shown in the discussion above is that this paradigm has a 
rigorous basis in quantum field theory. In closing, we thank Prof. Ignatios 
Antoniadis for the support and kind hospitality of the CERN TH Unit while 
part of this work was completed. 
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